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While the properties of the Kondo model in equilibrium are very well understood, much 
less is known for Kondo systems out of equilibrium. We study the properties of a quantum 
dot in the Kondo regime, when a large bias voltage V and/or a large magnetic field B is 
applied. Using the perturbative renormalization group generalized to stationary nonequilibrium 
situations, we calculate renormalized couplings, keeping their important energy dependence. 
We show that in a magnetic field the spin occupation of the quantum dot is non-thermal, being 
controlled by V and B in a complex way to be calculated by solving a quantum Boltzmann 
equation. We find that the well-known suppression of the Kondo effect at finite V ^ Tk (Kondo 
temperature) is caused by inelastic dephasing processes induced by the current through the dot. 
We calculate the corresponding decoherence rate, which serves to cut off the RG fiow usually 
well inside the perturbative regime (with possible exceptions). As a consequence, the differential 
conductance, the local magnetization, the spin relaxation rates and the local spectral function 
may be calculated for large V, B ^ Tk in a controlled way. 

KEYWORDS: Kondo effect, quantum dots, non-equilibrium, decoherence, perturbative renormalization 
group 



1. Introduction 

After 40 years of research, the single-impurity Kondo 
model^ and its variants are certainly among the best 
understood strongly correlated systems. An amazingly 
broad spectrum of powerful theoretical methods^ have 
been developed to describe the physics of a magnetic im- 
purity in a metal, ranging from exact analytical solutions 
using the Bethe Ansatz, bosonization and conformal field 
theory, diagrammatic methods, flow equations, pertur- 
bative renormalization group to sophisticated numeri- 
cal techniques like the numerical renormalization group. 
Using these tools an almost complete understanding of 
thermodynamic, transport and spectral properties of the 
single-impurity Kondo model has been achieved.^ 

There is, however, one important field about which in 
comparison little is known: Kondo physics out of equilib- 
rium, a question of both experimental and theoretical im- 
portance. Non-equilibrium can, for example, be reached 
when a finite current is driven through a Kondo system 
(see below). 

The study of the Kondo effect out of equilibrium 
started both experimentally and theoretically in the late 
1960's when the tunneling between two metals through 
insulating barriers was investigated'^"* as a function of 
bias voltage. The observed logarithmic enhancement of 
the tunneling for low temperature and small voltages was 
attributed to exchange tunneling through magnetic im- 
purities within the barriers by Appelbaum^' and An- 
derson. The relevant impurities were argued^^ to reside 
close to one side of the junction, and the localized spin 
was therefore tacitly assumed to be in equilibrium with 
the metal closeby. The tunneling to the other metal is 
weak and the conductance in such a situation is deter- 



mined by the equilibrium density of states on the im- 
purity. -'^'^ Subsequent work by Ivezic"'^^ considered also 
impurities deep inside the junction, which were pointed 
out to require a full non-equilibrium treatment. A good 
review of these earlier works may be found in Ref.^^ 

With the progress in nanotechnology, it became pos- 
sible to realize Kondo physics in quantum dots.^^"^^ A 
dot in the Coulomb blockade regime which carries a net 
spin can be mapped to the Kondo model. ^^'^^ The res- 
onant tunneling through the dot leads to a removal of 
the Coulomb blockade, i.e. an increase of the conduc- 
tance from small values up to the quantum limit, as the 
Kondo resonance develops. A finite bias voltage V drives 
the system out of equilibrium. It has been observed that 
the Kondo effect is quenched by raising the transport 
bias voltage V well above the Kondo temperature Tk, i.e. 
V ^ Tk , and that the presence of a magnetic field splits 
the zero-bias conductance peak into two distinct peaks, 
located at bias voltages roughly equal to plus and minus 
the Zeeman splitting of the spin on the dot.^^"^° Similar 
experiments were also possible in metallic nanoconstric- 
tions where it is possible to measure transport through a 
single magnetic impurity, however, in contrast to quan- 
tum dots it is not possible to control system parameters 
in such devices. 

The most straightforward way to study quantum dots 
out of equilibrium is to apply a finite dc bias voltage 
and to measure the current- voltage characteristics. There 
have, however, been a few remarkable experiments which 
go further. For example, Franceschi et al?^ managed to 
measure the splitting of the Kondo resonance by a dc bias 
voltage by using a three-lead configuration. Another set 
of questions can be addressed by studying the response to 
time-dependent external fields induced, for example, by 
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external irradiation by a microwave field with frequency 
^ 26,28-31 K^Qgan et al.^^ recently succeeded in observing 
satellites of the Kondo effect separated by ui in such an 
experiment. 

Kondo impurities also have a pronounced effect on the 
distribution of electrons in mesoscopic wires, which are 
driven out of equilibrium by an applied bias voltage. ^'^"^^ 
The inelastic scattering from the magnetic impurities at 
finite bias leads to a characteristic broadening^^"^^ of the 
electronic distribution functions which can be measured 
in tunneling experiments. ^^"^^ 

Theoretically, the Kondo effect out of equilibrium has 
been studied by a number of methods ranging from per- 
turbation theory,^' ^^^^ equations of motions 
and self-consistent diagrammatic methods^"' ^''"55 (using 
the so-called non-crossing approximation), slave-boson 
mean-field theories, ^''"^^ exact solutions for some vari- 
ants of the Kondo model with appropriately chosen cou- 
pling constants."''^ the construction of approximate scat- 
tering states starting from Bethe ansatz solutions, to 
perturbative renormalization group^^'^"^'^^ (reviewed be- 
low). It is, however, important to note, that many of the 
methods which have been so successful in equilibrium 
cannot or have not yet been generalized even to the sim- 
plest steady-state non-equilibrium situation. One of the 
reasons for this is that the current-carrying state at finite 
bias is a highly excited many-body state of the system. 
Therefore, all methods which by construction focus on 
ground state properties are not readily generalized to 
such a situation. 

In the following sections, we will review our version'*^ 
of perturbative renormalization group (RG) in the pres- 
ence of a finite bias voltage - formulated in the spirit of 
Anderson's poor man's scaling. We generalize the RG 
equations presented in Ref.^^ for a symmetric dot to an 
arbitrary exchange coupling of the spin to the conduc- 
tion electrons of two attached leads with different electro- 
chemical potential. The central goal is to perform a con- 
trolled calculation in the limit of either large bias voltage 
V ^ Tk or large magnetic field B ^ Tk or large probing 
frequency ui ^ Tk- We suggest a method to calculate the 
behavior of the conductance and other physical quanti- 
ties like magnetization or spectral fiuiction in leading or- 
der of the small parameter l/ln[max(y,i?,a;)/Tx]. We 
will argue that the perturbative renormalization group 
differs in three main aspects from its coimterpart in equi- 
librium:®^ First, the magnetization of the Kondo impu- 
rity (or in general the occupation probabilities of the 
quantum states) has to be calculated self-consistently 
from appropriate quantum Boltzmann equations.^^''*^''*® 
This leads to an unusual dependence e.g. of the spin sus- 
ceptibility on V and to a novel structure of the logarith- 
mic corrections. Second, the perturbative renormaliza- 
tion group has to be formulated in terms of frequency- 
dependent coupling functions instead of coupling con- 
stants. The reason is that electrons in an energy window 
set approximately by the external bias voltage contribute 
to the low-energy properties (sec Fig. 1) and their posi- 
tion within this window will affect their effective cou- 
pling to the spin. Therefore, the effective renormalized 
coupling of the conduction electrons to the local spin 



will depend explicitly on their energy. Third, decoher- 
ence effects are much more important out of equilib- 
rium.^^' ^^'^^'^^'^^ A finite current will induce noise and 
thereby induce dephasing of the coherent spin-fiip pro- 
cesses responsible for the Kondo effect. We will show that 
it is precisely due to those dephasing effects that a con- 
trolled calculation at large voltages is possible. 

2. Perturbation Theory 

We consider the Kondo Hamiltonian 

H = ^ (eic - A<a)c]jkCTCaka " BS^ 

+ ^ Ja'aS ■ icl,j^,^,T^,^Cak„), (1) 

a,a'=Z/,i?,k,k',(T,cr' 

where ^il.b. = i-V/2. S is the spin 1/2 on the dot and 
T are the Pauli matrices. This model describes a quan- 
tum dot coupled to two leads (the left and right electrons 
described by c^/^ j^^) in which the number of electrons 
on the dot is fixed to an odd integer by an interplay 
of gate voltage and charging energy. In this Coulomb 
blockade regime, effectively a single spin is localized on 
the dot which interacts via an exchange coupling J^a' 
arising from tunneling processes from and to the leads 
involving virtual excitation of the dot. In cases where 
the quantum dot can be described by a simple Anderson 
model, one obtains Jlr = Jrl = VJllJrr- In more 
complex situations, Jaa' can be an arbitrary symmet- 
ric 2 x 2 matrix. In general, a further cotunneling term 
Vcl^Ca' exists. As such a term does not fiow to strong 
coiipling under renormalization group, it can be safely 
neglected (furthermore, such a term vanishes for an ap- 
propriate choice of a gate voltage in the middle of the 
Coulomb blockade valley). We shall use the dimension- 
less coupling constants Qaa' = -^o-^aa', where Nq is the 
local density of states in the leads. For simplicity, we 
have coupled in Eq. (1) the magnetic field (measured in 
units where gfiB = 1) only to the local spin. An extra 
coupling to the electrons would effectively lead only to a 
small renormalization of the g factor. We assume a con- 
stant density of states in the leads in a frequency range of 
order B and V. In such a situation the density of states 
at the Fermi level remains unmodified in the presence of 
B. 

In order to be able to use standard diagrammatic tech- 
niques, we represent the local spin S = ^ '^^y f^'^'y-y'Iy 
by pseudo-fermions in the sector of the Hilbert space 
with /jj/^ = 1. Details of the Keldysh diagrammatic 
method®^ used throughout this paper and how the exact 
projection to the physical Hilbert space of the pseudo 
Fermions is performed can be found in Ref.^^ where the 
perturbation theory for the Kondo model is derived in 
detail to leading logarithmic order. 

Before setting up a renormalization group scheme, it 
is useful to investigate the results of perturbation the- 
ory, e.g. for the spin susceptibility. Out of equilibrium, 
the calculation of the susceptibility has to be done with 
some care, the bare perturbation theory diverges even at 
finite temperature T. Physicswise this arises because for 
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J — the magnetization of the completely uncoupled dot 
is undetermined. While in equilibrium any infinitesimal 
coupling leads to the usual thermal occupation indepen- 
dent of all details of the coupling, this is not the case 
out of equilibrium, where the details of the coupling do 
matter. Practically this implies that one has to calculate 
the occupation function even in the limit of vanishing 
J from some type of (quantum-) Boltzmann equation. 
In diagrammatics, this is given by one component of a 
(self-consistent) Dyson equation in Keldysh space (see 
Refs.^'*'"'^ for details). To lowest order this is just the 
well-known Boltzmann equation with golden-rule transi- 
tion rates. To obtain the steady state, the spin-flip rate 
from t to I has to equal the rate in the opposite direction 
and therefore 
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where f^j is the Fermi function. Solving for n| 
the limit B ^0. one finds''3.44,46 



X = 



91m 



(3) 



In equilibrium, y = 0, all coupling constants cancel out 
and one obtains x = 1/(2^) to zcroth order in the cou- 
pling constants. For a finite voltage, however, the result 
depends on the ratio of the coupling constants even for 
infinitesimal Jaa'- For large voltages one finds x ~ 1/^- 
One order higher, the logarithmic correction charac- 
teristic for the Kondo effect^ arises and one obtains, 
with ga = {gLL + 9rr)/'^, 



Fig. 1. Due to the finite bias voltage, the electrochemical poten- 
tial of the right and left lead differ by V . As all electrons in 
an energy window of width V contribute to low energy proper- 
ties, it is necessary to take into account that the effective ex- 
change coupling of the electrons depends on their energy. Right 
panel: frequency dependent effective coupling g{uj) calculated 
from Eq. (26) for V/Tk = 100. 



gle energy scale, the Kondo temperature Tk- At least 
for sufficiently small T/f, all microscopic complications, 
like band-structure effects, the energy dependence of the 
couplings, etc. can be absorbed in the value of Tk but 
will otherwise not affect low-energy properties which are 
completely universal functions of T/Tk or B/Tk for 
T,B ^ D. In technical terms, this arises because the 
Kondo model is renormalizable. A similar behavior is 
expected in the presence of a bias voltage V . We ex- 
pect e.g. that the conductance G through a quantum 
dot in the Kondo regime will be a universal function, 
G{V/Tk,B/Tk,T/Tk) at least for symmetric coupling 
Jll — Jlr — Jrr — J- (Note, however, that the con- 
ductance and all other physical observables will depend 



2<?: 



2 

LR. 



1 



25<ilnf 



2gd In 



Et- 577 ( 577 + 2^77 In f 



2giRln 



2V coth (^) glj, (l + In ^) + 2T g^^ (g^^ + 2g 



2 In^ 

77 T 



(4) 



where the bandwidth D cuts off the logarithmic singular- 
ities at high energies. Note that logarithmic corrections 
arise already to linear order in g with prefactors which 
again depend on the ratios of the coupling constants. In 
the limit of T 3> T^, in contrast, the logarithms all take 
the form ln(D/T'), the corrections in the numerator and 
in the denominator cancel, and we are left with the usual 
Curie law x = 1/2T', as in equilibrium the first logarith- 
mic correction arises only to order g^ , not included in 
Eq. (4). This shows again that in the presence of a fi- 
nite bias voltage, the structure of perturbation theory is 
strongly modified. 

3. Perturbative Renormalization Group 

Even for small coupling constants g and for sufficiently 
large V and B, such that y,i3 ^ Tr-, but still in the 
scaling regime V, B ^ D, such that g In D/{V,B) < 1, 
bare perturbation theory converges slowly. It is neces- 
sary to sum the leading logarithmic contributions in all 
orders of perturbation theory. Only after such a rcsum- 
mation another important property of the Kondo model 
becomes manifest: its physics depends only on a sin- 



on the ratios of the coupling constants for y > as is 
obvious from the simple fact that no current is fiowing 
for Jlr ^ 0.) 

Anderson^2 pioneered with his " poor man's scaling" a 
powerful method to resum the leading logarithmic cor- 
rections of the Kondo model in a controlled way: pertur- 
bative renormalization group (RG). This approach makes 
use of the fundamental idea that a small change of the 
cut-off D can be absorbed into a redefinition of the cou- 
pling constants g. As long as the cutoff-dependent run- 
ning coupling constant g{D) is small, the change oi g un- 
der an infinitesimal change of D, dg/dhiD, may be cal- 
culated in perturbation theory. In the equilibrium Kondo 
problem for vanishing magnetic field and temperature, 
the coupling constant grows when the cutoff is reduced, 
leaving the perturbative regime 5 < 1 when the running 
cutoff D is of the order of Tk ■ The ffow of the coupling 
constant to infinity leads finally to a complete screening 
of the localized spin. 

In the presence of a bias voltage V » Tk the situa- 
tion is more complex as we will show below: while the RG 
fiow is not cutoff by the voltage itself, it is stopped^'^' 
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Fig. 2. Diagrammatic form of the RG equation (6) . The strokes 
symbolize derivatives with respect to InZ? and a, a' = -L/R, 
(T, c' =T/J. and 7, 7' =t/i denote the quantum numbers of the in- 
coming and outgoing conduction electrons and pscudo fcrmions, 
respectively. Their frequencies are given by u)c, c<Jc! t^/j with 
Wc + a;/ = o;^ + a;^. 



before the strong couphng regime is reached by inelastic 

processes induced by the finite current through the sys- 
tem (at least for the experimentally most relevant case 

Jlr = Jll = Jrr)- 

To derive rcnormalization group equations in such a 
situation, one can for example start from a straightfor- 
ward generalization of functional renormalization group 
methods, described e.g. in Rcf..^''^ to Kcldysh diagrams. 
The basic idea is to track the evolution of all one-particle 
irreducible diagrams when an infrared cutoff is lowered. 
The advantage of this formulation is that it naturally 
includes not only frequency dependent vertices but also 
self-energy effects which are necessary to describe de- 
cohcrcncc. We will actually not follow this route here, 
but use a simpler approach more in the spirit of Ander- 
son's poor man's scaling by investigating more directly 
the scaling properties of diagrams. At each step wc will 
try to simplify the equations as much as possible: each 
approximation is tailored to be exact in leading order of 
the small expansion parameter 1/ ln[max(y, S, cij)/T/f]. 
A different and considerably more involved real-time RG 
scheme has been developed by Schoeller and Konig^^ but 
has not yet been applied to this problem. 

We start by analyzing the perturbation theory of var- 
ious physical quantities like the susceptibility (4) and 
ask the question whether a change in the cutoff D can 
be absorbed in a redefinition of the coupling constants. 
This turns out not to be possible: in Eq. (4) logarithmic 
corrections e.g. to gLR appear as 2hiD/V in the denom- 
inator, but as \n.D/V + \nD/T in the numerator. When 
the cutoff D gets smaller than V, different physical quan- 
tities (or numerator and denominator of the same phys- 
ical quantity) seem to require different renormalizations 
of coupling constants. This apparent contradiction to the 
principles of renormalization group is easily resolved by 
realizing that one should expect that the coupling con- 
stants depend on energy as explained above and in Fig. 1. 
For example, when analyzing the Boltzmann equation 
(2), one recognizes that the numerator in (3) arises from 
integrals which are confined to the vicinity of the Fermi 
energy, whereas in the denominator the energy integral 
covers a finite range of width V, and the couplings turn 
out to be quite different in the two cases [c.f. Fig. 1]. We 
therefore investigate in the following the origin of the 
frequency dependence of the vertex corrections. 

All leading logarithmic terms in perturbation theory 
stem from the simple vertex renormalizations shown in 
Fig. 2 when the real part ~ l/{u! ± B/2) of the pseudo 



fermion Green function is convoluted with the Kcldysh 
component of the electron line — 27riA^o tanh[((A; — 
/Xc«)/2r]. Using cutoffs symmetric with respect to Hl,r, 
respectively, one obtains at T = 
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where Alo depends on V, B and the incoming and out- 
going frequencies. Approximating the logarithmic deriva- 
tive by a Heaviside step function is valid in the two limits 
I Aw I ^ D and |Aa;| <C -D. To leading order in our small 
parameter 1/ ln[max(V, B)/TK\^ the detailed behavior at 
I Aw I D is not important. Similarly, it is sufficient to 
keep track only of the real parts of the coupling constants 

on one Keldysh contour denoted by g^'^^^yj^ for an 
incoming electron in lead a = L,R with energy Uc and 
spin cr =1, 1 interacting with a pseudo Fermion with spin 
7 and frequency w/ describing the local spin (see Fig. 2). 
Primed quantities refer to outgoing particles. A calcula- 
tion to higher order in l/ln[y/r/f] also has to take into 
account the imaginary parts and the full Keldysh struc- 
ture of the vertices which can be neglected here. 
From Fig. 2 and Eq. (5) we obtain 
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where 



Q^ = Q{D-\u:\ 



(6) 
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To simplify the notation, "*" represents in each term 
the one frequency, which is fixed by energy conservation. 
The initial conditions for Eqs. (6) at the bare cutoff are 



77'- 



The renormalization group equations (6) are rather 
complex and difficult to solve as each vertex depends on 
three frequencies. Fortunately, one can drastically sim- 
plify them to leading order in 1/ ln[max(V, B)/Tii\ using 
two approximations. First, the energy of the spin is well 
defined and the pseudo fermion spectral functions there- 
fore strongly peaked at w = ±5/2, which allows setting 
Uf (uj'f) to — 7-f ("I'^f )• Hence, one has to keep track 
only of a single frequency, the energy of the incoming 
electron. Second, the coupling functions appearing on 
the right-hand side of (6) depend only logarithmically 
on the frequency while the Heaviside step function is 
strongly frequency dependent. This allows to neglect the 
frequency dependence of the vertices by approximating 
/(Aw)e(i:>-|Aa;|) « /(0)e(D-|Aw|) on the right-hand 
side of (6). 

It is also useful to introduce a more compact nota- 
tion, which separates spin-flip processes g± from non-flip 
vertices g\\hy defining 
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The energy argument in the definition of g±_ has been 
chosen to simphfy the RG equations (sec below) . It turns 
out that the property gjj = is preserved under RG, 
therefore it is not necessary to keep track of a possible 
7 dependence of g^- Using these definitions and the ap- 
proximations described above one obtains 
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and for the longitudinal component 
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The initial values of the RG equations are (w) = 

The equations (9) and (10) together with the results 

for the dephasing rates and other physical quantities 
given below are the main result of this paper. In Ref.,^^ 
these equations have been presented for the special case 
gcia „ const, independent of the lead index, where they 
simplify drastically (see below). 

It is useful to discuss the symmetries of Eq. (9) and 
(10). Due to the hermiticity of the underlying problem 
one finds 



aa 
9±a 



(u;) = 
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(c). 



(11) 

(12) 



The coupling matrices simplify drastically in certain lim- 
its. For example, if the Kondo model is derived from an 
underlying Anderson model, then only a single channel 
couples to the dot and the matrix of bare couplings 
has only a single non- vanishing eigenvalue. This property 
is conserved under RG as the RG equations have the ma- 
trix structure dg = gMg and therefore one can express 
the matrices in LR space in terms of a single function 
9±/z{i^) 



g'l%{u)=2g^/,{co) 



COS( 

sin c 



(cos sin 0) (13) 



Due to (11), g± does not depend on a in this case, 

S-ltM =9±i{^) = a±{^)- 

In equilibrium, V = 0, it is well known that within one- 
loop RG, the two eigenvalues of the coupling constant 
matrix g = g"" do not mix. This is therefore also a 



property of the RG equations (9) and (10), which have 
in this case the matrix structure dg cx gg. For each of 
the eigenvalues of g one obtains for y = 

dgza{<^) 



dlnD 

dg±{Lj) 
dlnD 



= -2{g^{B/2)YQ^+„B (14) 
= -g^{B/2)g,{Q) + 60,-5/2] 



(15) 

In Ref.^^ we used these equations to calculate the spec- 
tral functions at large frequencies and magnetic fields. 
We checked that these equations indeed reproduce the 
perturbation theory for the T-matrix not only to order 
g(^ln[..] but also to order g"^ln^[..] as has to be expected 
from a theory which resums leading logarithmic correc- 
tions. Furthermore, we compared^^ our results to spec- 
tral functions calculated with numerical renormalization 
group and found that the relative errors are - again as 
expected - of order 1/ ln[max(u;, i3)/Tif]. 

For vanishing magnetic field B and spin-rotational in- 
variant bare couplings J, the renormalized coupling con- 
stants obviously do not depend on spin and the RG equa- 
tions simplify to 
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If the dot couples symmetrically to the left and right 
lead, one has 
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(17) 



If, furthermore, only a single channel couples symmetri- 
cally to the dot, one can use (13) to parameterize the 
renormalized vertices by just two functions g±{uj) = 
g±{—w) and fif^f (w) = 5zj.(— w) which are calculated from 
the RG equations 
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This limit has been discussed before in Ref."*^ 

How are physical quantities calculated from the renor- 
malized coupling constants? To leading order in our small 
parameter 1/ ln[max(V^, i?, a;)/rR-] it turns out that it is 
sufficient to replace in the 2nd order perturbation theory 
expressions (or, equivalently, in the golden rule expres- 
sions) the bare coupling constants by the appropriate 
renormalized couplings in the limit D — !■ 0. In each case 
we have checked this by comparing with perturbation 
theory to order g^, where the first logarithmic correc- 
tion arises.^"' Experimentally the most relevant quantity 
is the current 
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(19) 

where f^^^ = /(w =F ^^/2) are the Fermi functions in 
the left and right leads. The spin occupations n-^/i are 
calculated from the rate equation 



[see Eq. (2)] with n-f + nj^ = 1 and 
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A quantity of interest for spintronics applications is the 
spin-current'^^ Jspi„ = -^{Slz - Sr^), where Sl/r,z is 
the z component of the total spin of the conduction elec- 
trons in the left/right lead. 
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+ 4n^/«(l - /«_^^):^[5-_^(c^ - ^ 

Spin transport will be discussed in detail elsewhere.^* 

As we argued in the introduction, the RG flow will 
finally be cut off and controlled by the dephasing of co- 
herent spin-flips. Therefore one also has to calculate the 
dephasing rate F = I/T2 by replacing the bare couplings 
in the perturbative formula**^ by the renormalized cou- 
plings 



(23) 



7^ 
2 

This formula takes into account a partial cancelation of 
self-energy and vertex corrections. In a careful analysis of 
perturbation theory, it was shown in Ref.,^^ how decoher- 
ence rates cut off logarithmic singularities and therefore 
stop the renormalization group flow. It turns ouf*^ that 
different rates (combinations of 1/Ti and I/T2 in the lan- 
guage of NMR) enter the various logarithms. For the situ- 
ation discussed in this paper, however, the various relax- 
ation rates differ only by factors of order 1. As the relax- 
ation rates appear only (see below) in arguments of loga- 
rithms, InT/Tx, it is not necessary to keep track of these 
prefactors to leading order in 1/ ln[max(F, _B, ct;)/Tft:]. 
Hence, we can approximate all relaxation rates by F. To 
leading order, the effect that the self-energy and vertex 
corrections cut off all logarithmic contributions at the 
scale F, can be taken into account by replacing 



= e{D - ^/lo^ + r2) 



(24) 



in (9,10). A further effect of the relaxation rate is that 
spectral functions and spin-spin correlation functions are 



broadened on the scale F. This can effectively be de- 
scribed by replacing /^^ appearing in (19-23) by Fermi 
functions smeared out on the scale F (note, however, that 
the distribution functions of electrons in the leads are not 
renormalized). Details of how this broadening is justified 
and implemented will be published elsewhere. 

A full solution of the RG equation can now be obtained 
in the following way (an example is given below): First, 
the RG equations (9,10) are solved for F = 0. This turns 
out to be rather easy as on the right-hand side of the 
equations only couplings at a finite set of fixed frequen- 
cies enter. One therefore first solves for those frequencies 
and constructs the full energy dependence in the limit 
Z) ^ in a second step. For some frequencies one will 
find diverging couplings as the effect of F has not yet been 
taken into account. With these frequency dependent cou- 
plings, one calculates F from (23) which is finite as the 
integrals over the weakly divergent coupling functions do 
not diverge. With this F one recalculates the RG equa- 
tions using (24). In practical numerical implementations, 
we iterate this procedure until convergence is reached, 
but formally the effects of self-consistency are subleading 
and the first iteration described above is sufficient within 
the precision of our approach. Finally, physical quanti- 
ties like the magnetization [using (21)], the current (19) 
or the spectral function^^ are calculated from the renor- 
malized couplings. Note that the approximation, that F 
and other physical quantities like the occupations can 
be calculated in an independent second step after solv- 
ing the RG equations, is only valid to leading order in 
1 / In [max (y, B) /Tk] ■ To higher order, one has to solve 
simultaneously the RG equations for (frequency depen- 
dent) self-energies and vertices. 

As an example, we consider a simple situation by cal- 
culating the properties of a symmetric dot with g^^ = 
gRR _ fpj. vanishing magnetic field, B = Q, and 
large voltages, V 3> Tk- From (11) and (17) one finds 



(a;) 



g^^{uj) = g^^{-uj) = g'^{ijo) and g^^{uj) = g^^ 
g'^^{—Lo). To solve the RG equations (16), we first ob- 
serve that on their right hand side only the couplings 
gdiy/2) and g^^{V/2) enter. Therefore we first deter- 
mine the RG equation for these two constants 

-[g^'^(V/2fQv+g<i{V/2fQr] 



dgd{V/2) = 
dg'^\V/2) = 



5^^(y/2)ffd(i^/2)(ey + er) (25) 



with Or = e[Z)-F] and 6^ « Q[D-V] as F <C F. Note, 
that the RG flow is not completely cutoff by V but con- 
tinues for D > T. The equations are most conveniently 
solved in the scaling limit, where the Kondo temperature 
is held flxed assuming a large initial cutoff D ^ oo and 
a small bare exchange coupling J — » 0. For D > V, the 
RG equations are identical to the RG equations for the 
two-channel channel-anisotropic Kondo model in equi- 
librium, the channels are the even (L-|-R) and odd (L-R) 
linear combinations of electrons. The RG flow is char- 
acterized by two RG invariants, the Kondo tempera- 
ture Tk = De-i/tf^+f""") and ckq = ^^^r^. The RG 
equations are easily solved and one obtains gd{V/2) = 
for £) > y and gd{V/2) = 



l+2a„ ln[g/Tg] 



2ln[D/TK]{l+aoln[D/TK]) 
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l/\n[D/T*] with T* = T^^,(T^/v)i/(i+2aoin[y/rK]) for 
V < D < V . The couphng gLR{V/2) is given by 
gLR{V/2) = gd{V/2)/{l + 2ao ln[max(i^, V)/Tk]). Plug- 
ging these solutions into the right-hand side of (16), one 
can easily integrate the equations analytically for an ar- 
bitrary frequency. As the formulas are rather lengthy, 
we show only the result for ao = when only a single 
channel couples to the dot and 5^(0') = g^^{Lo) = 3(0;) 
according to Eq. (13). For D — > one finds 



5(a;)«^e[|a;-/3y/2|-y] 



+ @[V - Iw - l3V/2\] ^ 



A\n[\w- I3V/2\ITk] 



1 



ln[y max( 



PV/2\,V)/T, 
1 



where V 



V 



4in[y/T^; 

is calculated from (23) 



(26) 



In 

Fig. 2 a plot of g{Lo) is shown for V/Tk ~ 100, where 
r K, I.QTk- For V ^ Tk, the relaxation rate is larger 
than Tk and therefore one stays in the weak coupling 
regime, g{uj) <^ 1 for ao = 0. As described in de- 
tail in Rcf., for ao > 1 the perturbative RG breaks 
down for large voltages and one enters a strong coupling 
regime. (In that analysis'''''' we did not keep track of the 
frequency dependent couplings, nevertheless all results 
remain valid; g^n in Ref.^^ has to be identified with 

gLRi^ = 0)). 

Using (26) one can easily calculate physical quantities 
like the current. As the current is obtained from an in- 
tegral over g{uj) it is not sensitive to details of giuj) and 
one finds that the differential conductance is given by 
G Ri iehhi[v/TK]^ ' ^ result obtained before by Kaminski 
et al?^ A quantity like the susceptibility, however, is sen- 
sitive to the peaks in g{uj) and one finds (restoring the 
ao dependence) by solving the rate equation (20) 



X{V) = 



l + 2ao ln^(l-hao In^) 



V 



1 



(l-|-2aoln^)ln 



1" 7^(l+«o 1" T^) 



(1+ao lii^)ln^ 



(27) 



While in the presence of a finite bias voltage, large log- 
arithmic corrections arise [related to the peaks in (?(w)], 
one obtains in equilibrium to the same order of approx- 
imation, just X = 1/(27^)- Expanding this result in the 
bare couplings, one obtains Eq. (4). 

In Fig. 3, we show the renormalized coupling constants 
for various values of the running cutoff D in the presence 
of both a bias voltage V and a magnetic field B. While 
asymptotically the RG equations can be solved analyt- 
ically, the various crossover regimes have to be calcii- 
lated from a numerical solution following the approach 
described above. When the cutoff D is lowered, more and 
more scattering processes freeze out and only resonant 
contributions with uj « ±V/2 ± S/2 or w « ±V/2 ± B 
survive and lead to pronounced peaks in the renormal- 
ized coupling constants. 

Experimentally, the best accessible quantity is the dif- 
ferential conductance (j(V, B) as a function of bias volt- 
age and magnetic field. Our calculations can be applied 




Fig. 3. Renormalized coupling constants g^iio) (upper panel) 
and gz{uj) (lower panel) for B = IOOTk and V = 80Tk and 
for different values of the running cut-off D. A symmetric dot 
with bare couplings g"" = g is assumed. Peaks signal that res- 
onant scattering processes take place at those frequencies. The 
peaks are cutoff by the relaxation rate F. 



to all systems described by the simple Kondo model (1), 
e.g. quantum dots with an odd number of electrons in 
the Coulomb blockade regime, provided that either volt- 
ages or magnetic fields are large compared to the Kondo 
temperature, V, B 3> Tk, and provided that other exci- 
tations of the dot have a much higher energy. The lat- 
ter requirement is often not fulfilled in quantum dots or 
molecules. While a generalization of the above described 
methods to more complex dots and molecules with more 
levels is straightforward, this will introduce more fitting 
parameters. If the Kondo model (1) is valid, just two 
parameters are needed: the Kondo temperature Tk and 
the L/R asymmetry of the coupling Jll/ Jrr (assuming 
that only a single channel couples to the Kondo spin as it 
is the case in most single dot experiments) . Both param- 
eters can be determined by comparing the T dependence 
of the linear-response conductance G(T, V ^ Q, B — Q) 
to exact theoretical results determined from numerical 
renormalization group calculations.^" This allows for a 
parameter-free theoretical prediction of G{V/TktB/Tk) 
for large V and B. In Fig. 4 we compare''"^ our theo- 
retical results [determined from a numerical solution of 
Eqs. (9,10,19,21,23) using (13)] to experiments by Ralph 
and Buhrman,^'* where ratios of V/Tk and B/Tk larger 
than 100 have been probed. The peaks in the differen- 
tial conductance G{V/B,B/Tk) at V ~ B and their 
characteristic shape arises due to the interplay of several 
effects. For V ^ B and B » Tk, the spin is strongly po- 
larized, the Kondo effect is suppressed and transport is 
dominated by processes where the spin does not flip. For 
V > B a. new transport channel is opened, as the voltage 
provides a sufficient amount of energy to flip the spin. Si- 
multaneously, the magnetization on the dot changes and 
resonant spin-flip scattering from one Fermi surface to 
the other strongly renormalizes the effective couplings. 
The smooth drop at large voltages arises, because the 
Kondo effect is further suppressed by the voltage. The 
theory fits the experiment surprisingly well, especially 
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Fig. 4. Conductance measurements of Ref.^* (symbols) on metal- 
lic point contacts in magnetic fields 0.85, 1.7, 2.55 T (B = 
36, 72, 104T/f, with Tk « SOmK^*). Assuming that the corre- 
sponding point contact is described by a single-channel {Jlr = 
V JllJrr) Kendo model, ^ 4.2 is determined from G(V= 

0,S = 0,T = 50mK) = -^^JMLjRR^G^^^{T/TK), where G^y^ 

is known exactly from NRG calculations.^" This fixes all pa- 
rameters for our RG calculation (solid lines, T = 0). As the (B- 
dcpcndcnt) background is not known experimentally, wc subtract 
AG = Gb — 5.2 • W~^Goi^, where Gs is fitted to our results 
at large V . The experimental temperature T = 50mK leads to 
an extra small broadening at V = B. Note that the agreement 
between theory and experiment is best for large voltages and 
magnetic fields as expected. 



when one takes into account that a relative error of or- 
der 1 / ln[niax(y, B)/Tk] has to be expected from our cal- 
culation. Either the next-order correction is accidentally 
small or some error in the determination of Tk (we took 
the value quoted in the experimental papcr^**) acciden- 
tally improves the fit. It is, however, worthwhile to stress, 
that deviations of theory and experiments are strongest 
in the regime of small voltages and magnetic fields as 
expected from an approach based on an expansion in 
l/ln[max(V,B)/TK]. 

4. Conclusions and Outlook 

In the last 40 years the Kondo effect has motivated 

many new theoretical developments. It can be expected 
that the Kondo model will also play an important role in 
developing techniques to understand and describe quan- 
tum systems o\it of cqiiilibrium in regimes where pcrtur- 
bative methods are not available. For example, it would 
be extremely useful to establish novel numerical methods 
which are able to calculate e.g. the current-voltage char- 
acteristic of a quantum-dot in the Kondo regime. Besides 
the obvious practical use of methods to predict and ex- 
plain transport properties of single-electron transistors, 
molecules and similar systems, we believe that it will be 
also a conceptual challenge to develop the language and 
tools to describe and classify strongly interacting quan- 
tum systems in steady-state non-equilibrium. 

In this paper, we have discussed a simple approach to 
generalize perturbative renormalization group to a situ- 
ation, where a finite bias voltage is applied to a quantum 
dot in the Kondo regime. We derived the one-loop RG 
equations for a spin coupled to the electrons in two leads 
by an arbitrary matrix of exchange couplings. This al- 



lows a controlled calculation of spectral functions, differ- 
ential conductances, susceptibilities and other properties 
of the dot for large bias voltages and/or large magnetic 
fields and frequencies. The expansion is controlled by the 
small parameter l/ln[max(y, B,a;)/TR:]. For the future, 
it will be interesting to generalize this approach to more 
complex situations and to investigate systems where the 
perturbative renormalization group breaks down at large 
bias voltages. 
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